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Abstract 

We consider many-body effects on particle scattering in one, two and three dimensional Bose 
gases. We show that at T = these effects can be modelled by the simpler two-body T-matrix eval- 
uated off the energy shell. This is important in ID and 2D because the two-body T-matrix vanishes 
at zero energy and so mean-field effects on particle energies must be taken into account to obtain a 
self-consistent treatment of low energy collisions. Using the off-shell two-body T-matrix we obtain 
the energy and density dependence of the effective interaction in ID and 2D and the appropriate 
Gross-Pitaevskii equations for these dimensions. Our results provide an alternative derivation of 
those of Kolomeisky et al. |l|, || . We present numerical solutions of the Gross-Pitaevskii equation 
for a 2D condensate of hard-sphere bosons in a trap. We find that the interaction strength is much 
greater in 2D than for a 3D gas with the same hard-sphere radius. The Thomas- Fermi regime is 
therefore approached at lower condensate populations and the energy required to create vortices 
is lowered compared to the 3D case. 

PACS numbers: 03.65.Nk, 03.75.Fi 
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I. INTRODUCTION 



Recent experiments on the quasicondensation of a two dimensional gas of atomic hydro- 
gen |3| and the possibilities of confining dilute atomic gases in "low-dimensional" traps j|, |5| 
have stimulated interest in the possibilities of Bose-Einstein condensation in two-dimensional 
systems. It has long been known that, in the thermodynamic limit, Bose condensation is 
not possible in two dimensional homogeneous systems at any finite temperature because 
long wavelength fluctuations destroy long range coherence ||. Instead such a system un- 
dergoes a Kosterlitz-Thouless transition [|7| and acquires local coherence properties over a 
length scale dependent on the temperature — a "quasicondensate" M. In the limit T — > 
global coherence is achieved in homogeneous 2D systems and a true condensate then exists. 
In a trapped 2D system the modifications of the density of states caused by the confining 
potential enable a true condensate to exist even at finite temperatures ||. 

In most treatments of a Bose condensed gas in 3D, particle interactions are described 
by a delta-function contact potential whose strength is determined by the zero energy and 
momentum limit of the two-body T-matrix (T 2 b) which describes scattering in a vacuum. 
This leads to the standard form of the interaction potential (^A^f^a^/rnj <5(r), where a 3( j 
is the s-wave scattering length. At higher order it can be shown that the interactions are 



actually described by a many-body T-matrix (Tmb) 0) n, [12| which accounts for the fact 
that collisions occur in the presence of the condensate rather than in free space. In 2D 
this correction is critical because the 2D two-body T-matrix vanishes in the zero energy 
limit H, O, and thus we must include this correction (at least partially) even at leading 



order ||14|| . In this paper we develop an expression for the many-body T-matrix in terms of 
the two-body T-matrix evaluated at a shifted effective interaction energy. In one and two 
dimensions we obtain an effective interaction which depends on the energy of the collision, 
in contrast with three dimensional gases. 

The energy dependence of the effective interaction can be written as a density dependence, 
in which form the results can be applied to trapped gases. This leads to a Gross-Pitaevskii 
equation (GPE) describing the condensate wave function which no longer has a cubic non- 
linearity in ip, but instead goes as {ipl 4 ^ m ID and as (\ip\ 2 / In (V'l 2 )'^ i n 2D. Such a modified 
GPE has already been introduced by Kolomeisky ||T], [2|] and Tanatar ||15|| , using arguments 
based either on the renormalisation group or a Kohn-Sham density functional approach fL6|. 
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Our discussion in this paper is to show how essentially the same results can be obtained 
by a consideration of many-body effects on particle scattering and to relate this to well- 
understood treatments of the 3D Bose gas. Indeed, substantially the same treatment as used 
in 3D applied to the ID and 2D gases leads to the energy dependent effective interactions. 
The principle difference is that these effects must be taken into account in leading order, 
whereas in 3D they can be neglected in the simplest treatments and only become important 
at finite temperature or high density. 

In the following section we discuss the Gross-Pitaevskii equation, and the limits in which 



a system may be considered two dimensional. In section |TJ we then derive the many-body 
effective interaction for low dimensional gases, before considering its implications for ID 
gases in section |IV|. Finally, using this effective interaction we obtain a form of the two 
dimensional Gross-Pitaevskii equation, and we present the results of numerical solutions for 



both ground and vortex states in section VI. 



II. THE GROSS-PITAEVSKII EQUATION IN 2D AND QUASI-2D 

The macroscopic wave function for a Bose-Einstein condensate is found in mean field 
theory using a non-linear Schrodinger equation known as the Gross-Pitaevskii equation - 
where the non-linear term arises from interactions between the atoms of the condensate. 
Obtaining the form of the effective interaction in 2D, and describing its effect on the solutions 
of this equation are the main concerns of this paper. 

Currently, most BEC experiments have created three dimensional condensates, which are 
described by a GPE of the form 

- ^ V ^( r ) + WrMr) + N g 3d |^(r)| 2 ^(r) = ^(r), (I) 

where V tTa , p (r) is the external trapping potential, N Q is the condensate population, /i the 
chemical potential, and g3d is the coupling parameter describing the effective interactions. 
The coupling parameter is generally taken to be the zero energy and momentum limit of 
the two body T-matrix which in 3D is a non-zero constant g%d = Aittfasd/m, where a is 
the s-wave scattering length. The T-matrix has the contact potential form T 2 b(r,r') = 
g3d<5(r — r')<5(r) in the limit that all the momenta involved in typical collisions are much less 
than 1/R e , where R e is the range of the actual interatomic potential (which is not in general 
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equal to the scattering length a). 

The obvious extension of these experiments in order to achieve the goal of two di- 
mensional condensates is to confine a gas in an anisotropic trap such that the gas is 
tightly confined in the z direction. For a harmonic potential such a trap has the form 



Vtra P ( r ) = rnto 2 /2 (p 2 + z 2 /7), with l z = ^Jh/ 2mu z as the characteristic trap length in the 
tightly confined direction, where cu z = cj/7 1 / 2 . On decreasing l z (decreasing 7) the system 
will pass from being three dimensional to being two dimensional in a variety of senses. 

The system can first be called two dimensional once l z has merely been decreased suf- 
ficiently that the mean-field energy of the condensate is small compared to %u z . In this 
case the dynamics of the system in the z dimension are restricted to zero point oscillations. 
Nonetheless, if l z is still much greater than 03d, then two body collisions are hardly affected, 
and hence interactions can still be described by the three dimensional contact potential g^. 
Therefore, although in this case the third dimension can be factored out of the dynamics of 
the system, at short length scales the interactions are still three dimensional. This regime 
can be described using the 3D GPE of Eq. flU) with the assumption that the wave function 
can be factorised as 

, / x fmuj z \ l / A ( muj z 9 \ 

»<A*) = *M(-sf) exp (-^r 2 )- (2) 

Substituting into the 3D GPE, and integrating over z leads to a two dimensional equation 



^>{ P ) + ^ 2 /rV(p) + g'N |^(p)| 2 ${p) = p>(p), (3) 



where p = {x,y}, p' = p — hu z /2, and the coupling parameter g' is given by g' 3d = 
{moj z /2'Kfi) 1 / 2 g^. The subscript here refers to the three dimensional nature of the inter- 
actions whilst the prime indicates that g' 3d is a two dimensional quantity. 

The above factorisation of the wave function remains valid as l z is decreased further, but 
the assumption that the scattering is unaffected begins to break down when l z is not much 
greater than a 3c j. The effect of the confinement on particle interactions has been discussed 
in detail by Petrov and Shlyapnikov [17], [18], who found that a 2D contact potential can still 
be used but that the strength of the interaction becomes dependent upon the confinement. 
The coupling parameter which they obtained in this "quasi-2D" regime is 

^ 1/2 r 1 /muA 1 / 2 . {Bhu, 



5q2d 



m 



+ -^Z± In I 

a 3( j V 2-nn ) \ 2p7r 



(4) 
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where B « 0.915. This expression is valid when the conditions mg' 2d /2%h 2 , R e /l z , 2/j l /fiu 2 <C 
1 are satisfied. In the large l z limit the l/a^d term dominates and the scattering is three 
dimensional as considered above. However in the fully 2D limit the logarithmic term in 
Eq. OD) dominates and g becomes dependent upon \i. Equation was derived from solving 
the two body scattering problem within the potential causing the tight z confinement. We 
will now show how essentially the same result can be obtained in the fully 2D limit by a 
consideration of the many-body effects on particle scattering. 



III. THE T-MATRIX IN THE GPE 



In order to describe the interactions within a truly 2D BEC we must consider 2D scatter- 
ing in the presence of a condensate. This is described by a many body T- matrix Tub, and 
the coupling parameter which appears in the GPE is in fact given by the matrix element 
(k / |T]y[B(£')|k) evaluated in the limit of zero momentum and energy (k, k', K, E = 0). Note 
that the many-body T-matrix is in principle also a function of the centre of mass momentum 
K, but this will not be explicitly indicated in this paper for notational simplicity. This will 
not be important for the results presented since we will always take the limit K = in this 
paper. 

Before discussing the many-body T-matrix, however, we will first consider the simpler 
two-body T-matrix which describes collisions between two particles in a vacuum and for 
which analytical expressions exist . We will then show how the many-body T-matrix 
can be obtained from the two-body version in the limit appropriate for the study of BEC. 



A. The two-body T-matrix 



The two-body T-matrix describing scattering from an interparticle potential V^(r) is the 



solution to the Lippmann-Schwinger equation [20 



(k'|T 2b (£)|k> = (k'|V( ri -r 2 )|k) 

+ E(k'|V(r 1 -r 2 )|q) 



1 



E — (e 



sp 



+ e 



sp 



(q|T 2b (£)|k>, (5) 



'K/2+q T c K/2-q7 

where k and k' are the relative momenta of the two particles before and after the collision 
respectively, and K is the centre-of-mass momentum. The energy of a single particle state 



5 



is e S fF, where in the homogeneous limit e s ^ = h 2 k 2 /2m. The total energy of the collision 
is E and includes a contribution from the centre-of-mass momentum K which cancels the 
corresponding contribution from the single particle energies. The two-body T-matrix is 
therefore independent of K, as it must be in free space. 

The scattering event described here could be a single interaction (k'|V|k), or alternatively 
the particles may first make a transition to an intermediate state |q) (weighted by an energy 
dependent denominator) before interacting again to emerge in state |k'). The recursive 
nature of Eq. ([|) sums all possible processes for which |k) — > |k'). For many applications 
we only need the "on-shell" T-matrix where both the energy and momentum conservation 
laws are fulfilled. However, it is also useful to consider the more general off-shell form shown 
above, where the momenta and energy may take arbitrary values. 

It can be shown that, for interaction potentials of a finite range R e , the T-matrix is 
independent of the incoming and outgoing momenta (in the limit kR e ,k'R e <C 1) ||19|| . In 
the position representation this corresponds to an effective interaction which is proportional 
to 5(ri — r 2 ). This contact potential approximation is of great utility in solving the GPE 
where the zero momentum limit of the T-matrix is used to describe particle interactions. 
In the three dimensional case the T-matrix elements at low energy and momenta are also 
independent of energy, leading to a constant coupling parameter in the GPE with form 
g 3d = 4:TTh 2 a 3d /m in Eq. ([[]). 

The contact potential approximation is still valid in one and two dimensions, but the 
T-matrix at leading order now depends upon the energy of the collision, as will be shown 
in the following sections. Thus the scattering terms in the 2D GPE will be quite different 
from the three dimensional case. 



B. The many-body T-matrix 

The two-body T-matrix describes collisions in vacuo in which the intermediate states are 
single-particle in nature. However, in a Bose condensed gas collisions occur in the presence 
of a condensate and a many-body T-matrix is needed to describe scattering processes. This 
is defined by the equation 

(k'|T MB (£)|k) = (k'|V(r 1 -r 2 )|k> 
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+ Uk'W^ - r 2 )|q> % + t " K ^-> (q|r MB (g)|k), (6) 

q — l^K/2+q + ^K/2-qJ 

where is the interaction energy, and e q is the energy of a quasiparticle state of momentum 
q, which is given by 

e 



'V 



4 P ) 2 + 2£ >f /2 > (7) 



in the Bogoliubov approximation |2T| for the case of the hard sphere gas. The corrections 
included in this many-body T-matrix over the two-body version are the occurrence of quasi- 
particle rather than particle energies for the intermediate states, and the Bose enhancement 
of scattering into these states. This latter effect results in the presence of population factors 
n q in Eq. (§). 

Formally, this many-body T-matrix is included in the theory of a Bose condensed gas by 
considering the effect of the so-called anomalous average (aiCij) on the condensate evolution, 
where hi is the non-condensate annihilation operator for state i. This term occurs when 
terms in the Hamiltonian of higher than quadratic order in Sj, a\ are taken into account [II], 
29|. We note that a generalisation of the many-body T-matrix which includes quasiparticle 
propagator factors for the intermediate states has been proposed JTIJ, but the corrections 
this includes over and above Eq. (§) are of still higher order. 

We note that the energies e q and E in Tmb are measured relative to the condensate, 
whereas the single particle energies in T 2 h are measured relative to the energy of a stationary 
particle. This means that for collisions between particles in the condensate we take the limit 
E = in Tmb; which corresponds to E — 2{i when measured relative to the same zero of 



energy as the two body case |10) . For collisions between condensate atoms, we also take the 
zero momentum limit k, k', K = 0. Interactions between two condensate atoms are therefore 
described by the matrix element (0|Tmb(0)|0). 



C Tmb in terms of T 2 h, a simple argument 

The Lippmann-Schwinger equation for the many-body T-matrix is substantially more 
difficult to solve than the two-body equivalent due to the presence of quasiparticle energies 
and populations. In the limit of zero temperature we will show that the many body T- 
matrix can be approximated by an off shell two body T-matrix evaluated at a negative 
energy. To see this we consider Eq. for the matrix element (0|Tmb(0)|0) at T = 
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where the population terms vanish. Upon comparison with Eq. (||) it can be seen that the 
only difference between the equations for the two types of T-matrix occurs in the energy 
denominators. Specifically, the quasiparticle energy spectrum appears in the many-body 
case, whereas the single particle spectrum appears in the two-body case. Heuristically, if 
the dominant contribution to the intermediate states in a collision comes from states with 
energies of order \i or higher, we can proceed by replacing Sk by e s k p + /i. This is the high 
energy limit of the Bogoliubov spectrum of Eq. (0) and it contains a constant shift from 
the single particle spectrum due to the mean field effects of the condensate which do not 
vanish in the relevant momentum range k ~ k$ for a contact potential interaction (where 
h 2 k 2 /2m = //). We are interested in the many-body T-matrix at E = 0, and thus the energy 
denominator in Eq. becomes 

1 _ 1 

- (>K/2+q + EK/2-q) - (e£ /2+q +H + ^K /2 _ q + V) 

— 2[1 — (^K/2+q + £ K/2-q) 

Comparison with Eq. (||) for the two-body T-matrix shows that in this approximation 

(0|T MB (0)|0) = (0|T 2b (-2^)|0). (9) 

Interestingly, this shows that the effective two body interaction energy is negative, meaning 
that the interaction strength is always real. We will see that this is important in the ID 
case in the following section. In 3D the two-body T-matrix is independent of energy to 
first order, but in both one and two dimensions it has a non-trivial energy dependence and 
therefore the effective interaction energy becomes important in these lower dimensions. 

At first glance the result of Eq. @ may appear counterintuitive since the energy of a 
collision between two condensate particles might be thought to be +2/z, and certainly not 
negative. However, as we have shown, the many-body effects in the system lead to a shift 
in the quasiparticle energy spectrum and it is this which leads to a shifted effective energy 
entering the two-body T-matrix. Stoof et al. |22| , have also proposed that interactions 
in low-dimensional condensates can be described by the two-body T-matrix evaluated at 
a negative energy (—2fi), the same result given by our heuristic argument above. In the 
following section we will use a more rigorous argument and find that this leads to somewhat 
better values for the effective interaction energy. 



S 



D. Tmb in terms of T2t>, a better argument 



Having shown heuristically in the previous section that the many-body T-matrix can be 
approximated by a two-body T-matrix evaluated at a negative energy, we will now present 
a more formal justification. This will lead to a slight modification to the magnitude of 
the energy used in the two body T-matrix, but the essential physics of the argument is 
unchanged. 

From equations (^) and (^|), it is possible to derive an expression for the many-body 
T-matrix solely in terms of the two-body T-matrix [|TIJ 



(k'|T MB (£)|k> = (k'|T 2b (£)|k> + (k'|T corr (£,£)|k>, (10) 

where 

ru\rr (w PMLA V" (k / |r 2b (E)|q)(l + nK/2+q + n K /2- q )(q|TMB(^)|k) 
{k\l corr {E,E)\k) = ^ — ■ r 

q^O ^ ~~ l £ K/2+q + ^K/2-qJ 

v (k / |T 2b (E)|q)(q|r M B(g)|k) 

^ F-(r sp 4- p sp \ ' ^ ' 

q ^ l £ K/2+q + £ K/2-qJ 

If we now assume that there is a value of E = E* for which (k'|TMB(£')|k) = (k'|T 2b (i?)|k), 
we can replace Tmb(E) on the RHS of Eq. ([TT|) by T 2 h(E*). The value of E* may then be 
found by solving for (k'|T corr (£', E*)\k) = 0. We again take the limit of zero temperature, 
such that nK/2+ q) %/2-q are zero, and for collisions between two atoms in the condensate 
we take the limit k, k', K, E = 0. The value of E* is then given in D dimensions by the 
solution to 

= / dk - / _ 2 - - dk. (12) 

Jo -2e k Jo E*-h 2 k 2 /m V ; 

Substituting the Bogoliubov dispersion relationship for the quasiparticle energies using 
Eq. (^) and carrying out the integrals in Eq. (12) we can obtain expressions for E*. We 



are then able to express the coupling parameter which occurs in the GPE in terms of the 
two-body T-matrix evaluated at the energies E*. In two and three dimensions this leads to 

f (0|T 2b (-u)|0) in 2D, 
9=0^0 0= ^ (13) 

[ (0lT 2b (-i| M )|0> in 3D. 

However, in ID the situation is more complicated because the first integral is logarithmically 
divergent. This case will be dealt with in section [IV| where we show that the results obtained 
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are consistent with known exact results. Prior to that however, we derive in the following 
section the form of the many-body T-matrix in two dimensions and show that the effective 
interaction energy becomes important in this case. 



E. The many-body T-matrix in two dimensions 

We consider the case of a 2D Bose gas with an interatomic potential V(p) which is short- 
range, parameterised by a length a 2d , arid which admits no bound states. Specifically, we 
consider the case of a 'hard-disc' potential such that V(p) = oo for \p\ < a 2d and V(p) = 



otherwise. In recent work we have derived a full expression for the two-body T-matrix 



for this potential in the general off shell case. In the limit ka 2d , k'a 2d "C 1 the result is 

<k'|T 26 (£)|k> = . - 4? f^ 2 7 - a r , (14) 
7rz — Z7 EM — m(Ema2 d /8n ) 

where the corrections are of order (ka 2d ) 2 and 1/ [ln(/ca 2 d)] 3 or greater, and 7em is the 



Euler-Mascheroni constant. This result agrees with the work of Stoof [24], and also in the 



half-on-shell limit with the earlier work of Schick |Tjj and Bloom ||25|| . It is also of the 
same form as the results obtained by Fisher and Hohenberg [|X3j] who considered the case of 
a Gaussian interatomic potential, implying that the result is general for most short-range 
repulsive potentials which may be parameterised by a length a 2 d- in this low momentum 
limit the T-matrix is independent of both k and k' and thus it is still represented in position 
space by a delta function effective interaction potential. The new feature compared to the 3D 
case is that the T-matrix now depends on energy and in particular it vanishes as E — > 0. It 
is therefore crucial to take into account the many-body shift in the effective collision energy 
of two condensate atoms. This is now a self-consistent problem as many-body effects give 
rise to a non-zero coupling constant. In 3D the two-body T-matrix is non-zero as E — ► 
and many-body effects can therefore be neglected at leading order for dilute gases. 



From equations (|H|) and (14) the many-body T-matrix, and therefore the coupling pa- 



rameter, in 2D is found to be 

, , m / \ i \ Airh 2 1 
S2d = T 2d (-p = — — — Y , 15 

where terms of order 1 / [ln(pma 2d / Ah 2 )] 2 or greater have been neglected. Note that the eval- 
uation of the two-body T-matrix at a negative energy means that the imaginary component 
in Eq. fll4|) vanishes, and thus the many-body T-matrix is real. 

10 



The parameter which appears in this description of the interparticle interactions is the 
two dimensional scattering length a 2( j, analogous to the 3D s-wave scattering length a 3( j 
which parameterises three dimensional collisions in cold dilute gases. Reliable values of 
have been obtained in 3D by experimental measurements, and potentially a 2 d could be 
measured in this manner. However, in their work on quasi-2D scattering processes Petrov 



and Shlyapnikov also derived an expression for this parameter in terms of the three 
dimensional 03d, and the confinement of the trap in the tight direction l z 

«2d = 4,/^/ z exp ( -y/n— ) , (16) 
V B \ a 3d J 

where B « 0.915. Using this expression for a 2 d in Eq. (|15|) we obtain the quasi-2D coupling 
parameter of Petrov and Shlyapnikov given in Eq. (£|), and our approach therefore agrees 
with their results in the genuine 2D limit which is appropriate for l z < ca- 
using this expression we are able to compare the strength of the 2D and quasi-2D cou- 
pling parameters with the parameter for quasi-2D gases with 3D scattering g' 3d described in 
section [TJ . These quantities are displayed as a function of trap width in the z dimension 
in Fig. |]. It can easily be seen that the size of coupling parameter appearing in the GPE 
for the genuine 2D case is over an order of magnitude greater than in the case where the 
scattering is essentially 3D in nature (l z /a>z& 3> 1)- The magnitude of g 2 d decreases slowly as 
l z is decreased beyond ~ (13d /2 (not shown on the graph) due to the size of a 2 d determined 
from Eq. (|16|), and it matches the 3D scattering limit for l z /a 3d > 10. 



IV. THE MANY-BODY T-MATRIX IN ONE DIMENSION 

Before we use the result in the previous section to solve the two dimensional GPE, we 
briefly consider the one dimensional case. Our discussion in this section is not intended to 
be rigorous, but is meant instead to demonstrate the importance of including many-body 
effects, via the many-body T-matrix, when considering the properties of a Bose gas in low 
dimensions. 

A one dimensional condensate is described by the Gross-Pitaevskii equation 

h 2 

- —V 2 ^{x) + V tTap (x)i>(x) + N o9ld |^(x)| 2 ^(:r) = ^(x), (17) 

where gid is the one dimensional coupling parameter. The use of the GPE necessarily 
assumes the existence of a condensate, which in one dimension implies that the system must 
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FIG. 1: Log- log graph of the effective 2D interparticle interaction strength as a function of con- 
finement in the third dimension. The solid line shows 2d which describes scattering in quasi-2D 



gases, taken from ref. [18]. Our results for <72d derived in this paper are consistent with this result 
and were derived for the region of validity shown. The dashed line shows g' sd which is the expected 
limit at large Iz/azd- 



be confined in a trap and therefore of a finite size. In a homogeneous ID system a true 
condensate may not exist in the thermodynamic limit due to the density of states || 



With this caveat in mind we will use the ID case to illustrate the importance of the energy 
dependence of the many-body T-matrix. Specifically we will consider the one dimensional 
analogue of a hard sphere gas for which exact results exist. This gas has an interatomic 
potential of the form 

for \x\ > 
oo for \x\ < aid- 



V(x) 



(18) 



In a recent paper [19] we have used an inhomogeneous Schrodinger equation to obtain 
results for the general off-shell two-body T-matrix for hard sphere gases in one, two, and 
three dimensions. In one dimension in the limit of zero momenta the result is 



(0|T 2b (£)|0> 




iy/Eaid + Ea 2 ld ) for E > 



E\a ld -Eai d ) for E < 0. 



(19) 



As in the two dimensional case, the T-matrix is dependent on the collision energy even at 
lowest order, and so the shift to an effective interaction energy predicted in section [111 C 



due to many-body effects is again important. Furthermore, in the case that Ea\ d <C 1, 
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the leading order term in the T-matrix in ID is imaginary if E is positive. The shift to a 
negative effective interaction energy is therefore critical in this one dimensional case. 

In order to obtain the many-body T-matrix in terms of this two-body T-matrix we must 
solve Eq. (|T2]). As noted earlier the first integral in this equation is logarithmically diver- 
gent in ID. Physically this arises from the fact that a true condensate does not exist in a 
homogeneous ID system. Instead of a single quantum level with a macroscopic occupation 
(as occurs in a true condensate), in ID there is a band of low energy levels which all have 
large occupations. The same methods as discussed above may still be used in the ID case 
however, provided that we now define the "condensate" as a band of levels in momentum 
space up to a cut-off at fc max , such that the ID "condensate density" n = Z)fe.<fc max n i sat- 
isfies n Q ~ n (as for a true condensate). Using this definition, the lower limit of the first 



integral in Eq. (12) should then be fc max and the divergence is removed. This approach is 
justified for a confined ID system since we may assume the existence of a condensate due 
to the modification to the density of states which also removes the divergence. 

A reasonable value for /c max may be obtained from the momentum distribution for a 



system of impenetrable bosons. Such a distribution is discussed in Ref. [2?|. We will define 
^max by the criterion that N{k > fc max ) < 1, which gives fc max ~ 0.257m pTfl . Using this as 
the cut-off in Eq. (0), a solution may be found for E* by making the ansatz that E* = —Cfi, 
where C is a constant. For the hard sphere case considered here the ansatz is satisfied when 
C is the solution to 



tanh 1 



TV 2 



(20) 



/ (0.25tt) 2 /4C + 4_ 

This can be solved numerically to give C ~ 3.4. The expression for the many-body T-matrix 
in ID to leading order is then 



Old = (0\T 2h (-Cfj,)\0) = ^ACh 2 ^/m. (21) 

We now consider a homogeneous ID Bose gas, using the above definition of the con- 
densate. When this system is in the ground state the contributions in the GPE from the 
curvature of the wave function and the trapping potential both vanish, and therefore 

ACh 2 2 . , 

A*id = n Q gi d = n . (22) 

m 

This form differs from that found in 3D where /i cx no because of the dependence of gid on 
the chemical potential. 

13 



This result can also be explained heuristically, as the extra curvature introduced into the 
wave function by the presence of the other atoms. If we consider a many-body wave function 
which scatters off a hard sphere potential of range aid, then in the limit of zero energy, we 
need to solve cPi/j/dx 2 = 0. We impose the boundary conditions that ip(x) = at x = aid 
and ip(x) approaches an asymptotic value x a t large x. Since ip is a many particle wave 
function, the distance at which it must arrive at its asymptotic value will be of the order of 
the interparticle spacing 1$. This gives a solution to the scattering problem of 

= 7 ( x ~ fl id) f° r a ld < x < l . (23) 

'0 — Old 

The extra energy caused by the curvature of the wave function in this region is then 

h 2 r l o h 2 M 2 

1 \V^j{x)\ 2 dx ^ -— (24) 



2m J a 2m/o 
And since Iq = 1/no and \x\ 2 = no we have that the interparticle interactions make a 
contribution to the energy which scales as rig. The same result may be derived from an even 
simpler argument which considers each particle to be confined in an infinite square well of 
length ~ 1/n by its nearest neighbours. 

The exact result for /i in such a ID gas has long been known. In solving the system of ID 
interacting bosons by demonstrating equivalence with a gas of ID non-interacting fermions, 



Girardeau [28| showed that in the strong coupling limit (appropriate to the hard-sphere 



potential considered above) 

fl =--n 2 . 25 
2 m 

Our result therefore shows the correct dependence on n 2 , but disagrees on the numerical 
factor. The disagreement is due to the fact that, as previously mentioned, in a homogeneous 
ID system there can never be a true Bose condensate, so significant corrections to the 
GPE can be expected. The additional uncertainty in the choice of fc max also introduces a 
source for discrepancy in the numerical factor. However, the agreement with the dependence 
on n 2 indicates that the energy dependent many-body T-matrix appears to deal with the 
interactions correctly. This is interesting because it means that an intrinsically many-body 
effect, namely particle confinement by neighbours, can be modelled by an off-shell two-body 
T-matrix evaluated at a shifted effective interaction energy which is the essential argument of 
this paper. This suggests that the method will have at least qualitatively the correct density 
dependence in the strong coupling limit. A more detailed investigation of this approach in 
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the ID case will be the subject of a further paper. Although our discussion in this section 
has been qualitative due to the lack of a true condensate in a ID homogeneous system even 
at zero temperature, in a trapped ID system it is possible for a true condensate to form. We 
therefore expect that semi- quantitative results outside the normal BEC regime of validity 
can be achieved using this method. In two dimensions a true condensate can be formed, 
even in a homogeneous system, at T = and so we expect our 2D results in this paper to 
be quantitatively correct. 



V. SCATTERING IN INHOMOGENEOUS GASES 



In the previous two sections we presented expressions for the many-body T-matrix in 
one and two dimensions in terms of the two-body T-matrix evaluated at shifted effective 
interaction energies. However, the results obtained are strictly only valid for homogeneous 
systems since we have not accounted for any modifications of the scattering wave functions 
due to the presence of a confining potential. We consider here the case of a gas confined 
tightly in one or two dimensions (in order to reduce the dimensionality, as discussed in 
section [II]) and weakly in the remaining dimensions on a length scale /trap- 

Provided that the range of the interatomic potential R e is much smaller than Z t rap then 
the scattering will be locally homogeneous and we can replace \l where it occurs in equa- 
tions fllBD and ( [21] ) by the homogeneous expression fi = nog. This is a form of local density 
approximation and, as the density of an inhomogeneous gas is spatially dependent, this 
leads to spatially dependent coupling parameters. Recognising that rao(r) = JV |^(r)| 2 the 
coupling parameters in one and two dimensions are 



9id = mx)\ 26 

m 

«- = ~[^WI^)]-' + <.(!^l) (-) 



m 

These results agree with the work of Kolomeisky el al. [|l], § who obtained similar expressions 
based on a renormalisation group analysis. Such density dependent coupling parameters are 



also expected from the results of density functional theory [16] which predict that the energy 
of the system is a functional of the density only. The same results may be obtained from 
mean-field theory by incorporating the spatially dependent anomalous average (ctiCij) into 



the system of equations governing a condensate and solving self-consistently |29|j . 
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VI. 2D SOLUTIONS OF THE NON-LINEAR SCHRODINGER EQUATION 



In this section we present solutions of the GPE for a trapped two dimensional gas. The 
solutions are found for a given p by propagating the time dependent GPE forward in imagi- 
nary time from an initial approximate solution to obtain both the ground state wave function 
and the non-linearity #2d N . As mentioned in the previous section the coupling parameter 
in two dimensions in a trap is spatially dependent, having a logarithmic dependency on the 
density. However, since in two dimensions the spatial dependence is merely logarithmic it 
will have little effect on the solutions of the GPE, except at the very edges of the trap where 
the wave function vanishes. We therefore use the homogeneous system coupling parameter 
of Eq. (|T5D, which will illustrate the features of most interest. 

Using the expression for the 2D coupling parameter found in Eq. ( |T5D we solve the two 
dimensional time-independent Gross-Pitaevskii equation for a 2D Bose condensate in a trap 
with V tTSL p(p) = ^mLJ 2 p 2 . We can make the GPE dimensionless, scaling all energies by Tiuo 
and all lengths by l p = h / (2mu) , giving 

- VV(P) + W/Mp) + Nofc d (p) mp)\ 2 ^(p) = Mp), (28) 

where #2d(p) = — 87r/ ln(//a| d /8) and Vt rap (p) = \p 2 . Note that the quantity pd 2 d is small 
compared to unity (or the earlier expansion of the T-matrix elements fails) and therefore the 
interaction is repulsive. As shown earlier, for the range of p and N that we consider here, 
we find that Eq. (|T5| ) leads to a value for <jbd(p) which is more than an order of magnitude 
greater than the equivalent value for a quasi-2D gas in which the particle interactions are 
effectively 3D in nature. Thus the non-linear term in the GPE is more significant in 2 
dimensions than in the 3D case. 

A. Ground state solutions 

Figure || presents sample solutions for the ground state of a 2D BEC in a trap for differing 
values of p. To illustrate the physical quantities involved we give numbers for a gas in a 
trap of uj = 2n x 100Hz and with a scattering parameter given by a 2 d = 6nm. This is 
close to the 3D s-wave scattering length found for 87 Rb, and therefore from Eq. (|16|) 
this corresponds to a situation where l z ~ a 3d . We see that at low N the solution is 
approximately the Gaussian wave function that is expected for the non-interacting case. At 
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FIG. 2: (a) Ground state 2D GPE solutions in an axisymmetric trap for \l = 2, 10, 20, 30, 40, and 
50Tkj. ip(p) is normalised to unity, and populations given assume parameters to = 2irx 100Hz, a>2d = 
6nm. (b) Comparison of GPE solution for \x = 50 (solid line) with Thomas- Fermi approximation 
(dashed line). 

higher N the Thomas-Fermi approximation found by neglecting the contribution to the 
GPE from the kinetic energy term as compared to the interaction and trapping terms is 
expected to be a good description. In two dimensions the Thomas-Fermi approximation 
gives a density profile in the form of an inverted parabola 



TF\ 



hl(/^2 d /8) 



iV 87T 



fi ~ Krap(p)J 0(ji- Krap(p)J , (29) 

where 8(x) is the step-function. At higher N the solutions shown are generally very well 
approximated by the Thomas- Fermi form, except at the boundary region of the condensate. 
Indeed we find that the Thomas-Fermi approximation works well for the 2D case, due to the 
high strength of the scattering, as expected from the dimensionless GPE fl28|). For the non- 
linear term to dominate the kinetic energy term requires that iV^ 3> — ln(/xa| d /8) (where 
p, ~ 10 — 100), whilst in the three dimensional case we require NQ d ^> l/a^. Putting typical 
numbers into this using our parameters we get N^ d ^> 10, whilst ^> 100, and thus the 
Thomas-Fermi regime is reached in 2D with about an order of magnitude fewer atoms than 
is the case for 3D. As confirmation of this, the Thomas-Fermi approximation for the number 
of condensate atoms is N = — jj 2 \n \ uma\ A / 'Ah 2 ^j /{2huo) 2 and is found to be in agreement 
with the numerical results to within one percent once N was greater than 300. 
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FIG. 3: Nq vs /x for a 2D Bose gas. The dots represent solutions of the GPE with the full energy 
dependent interaction given in Eq. (U). The lines are results which assume a constant (independent 
of /x) coupling parameter g 2 &. The three constant values of g 2 & correspond to Eq. ( |15D evaluated 
at \x equal to 6Hlo (dotted), 25hu (dashed) and 50huj (solid). 



In some previous papers |30| the GPE has been solved with g 2( i approximated by an 
energy independent constant. This is appropriate to the case where the scattering is three 
dimensional, but not to the fully 2D case where g 2 d depends on /x. We find here that the 
interaction strength given by Eq. flI5| ) increases by about 50% as /x rises from 2hu to 50hu. 
Figure |3] shows the possible errors which can arise from making the assumption of a constant 
coupling parameter. Each line plotted on this graph assumes a constant g 2 d, the strength 
of which is chosen to agree with Eq. ( |15|) at a certain value of the chemical potential /x*. 
The figure shows that results obtained with a constant <?2d (/•**) will introduce systematic 
errors when /x is significantly different from /x*. In the Thomas- Fermi approximation the 
relative error incurred in a measurement of iVo assuming a constant g 2 d{p*) is given by 
ln(/x*//x)/ln(/xai d /8). 

B. Vortex state solutions 

The 2D GPE can also be solved for the case of a two-dimensional condensate in a sym- 
metric trap containing a vortex at the centre by looking for solutions of the form 

ij(p) = <P(\p\)e iK6 , (30) 
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FIG. 4: Sample 2D GPE solutions for a vortex state with k = 1 and for values of p of 3, 5, 10, 15, 
and 20Huj. (f)(p) is normalised to unity, and populations given assume parameters u = 2ir x 100Hz, 
a 2d = 6nm. 

where 9 is the angle around the vortex core, and the phase wraps around by 27tk, where k 
is an integer, as the range of 9 is traversed. This adds an "effective potential" to the GPE 
and we now solve 



Wp)0(p) + N g 2 M |0(r)| 2 0(p) = p0(p). 



(31) 



2m ' ' ryr/ ' 2mr 2 

Solutions of these vortex states are shown in Fig. 

Such vortex states, which carry an angular momentum L z = Nhn, can be made ener- 
getically favourable by rotating the trap with sufficiently high frequency Q. The energy 
functional for a wave function in the non-rotating frame is 



E\ 



jd P 



|^|v^(p)| 2 + Wp) Wp)| 2 + ^%(p)I 4 



(32) 



The point at which E[ifj K= i] — QL Z becomes less than E[if) K=Q ] is known as the ther- 
modynamic critical frequency, and this is plotted in Fig. [| for both 2D and (genuine) 3D 
condensates. The three dimensional results were calculated from solutions of the 3D GPE, 
given by Eq. (|l|), with a 3 a taken to be equal to a>2d the scattering length used for the 2D 
results. Creation of a vortex in the centre of the trap comes at the cost of increasing the 
contributions from both the kinetic energy and the trapping potential terms in the GPE, 
although the non-linear contribution is reduced by virtue of a lower central density. Stronger 
nonlinear systems are therefore more susceptible to vortex creation, and this becomes ener- 



19 



1 



u 0.4 

a 




0.2 












2 



4 



6 

10- 3 iV( 



8 



10 



12 



o 



FIG. 5: The critical frequency O c at which vortex formation becomes energetically favourable in 
2D (lower) and 3D (upper) gases as a function of condensate population Nq. Results obtained 
using a 3d = a 2 d- 

getically favourable at much lower frequencies in 2D than in 3D for the same value of the 
scattering length a, as seen in Fig. |^. 

VII. DISCUSSION AND CONCLUSIONS 

In this paper we have found expressions for the many-body T-matrix in a dilute Bose 
gas describing the collisions occurring in a condensate in terms of the simpler two-body 
T-matrix. We have shown that many-body effects of the condensate mean-field on such 
collisions may be incorporated by a shift in the effective interaction energy of a two-body 
collision, and that such an approach leads to the same results obtained from renormalisation 
group techniques fj. 

The fundamental difference to the three dimensional case is that the first order term in 
the T-matrix in lower dimensions is dependent not only on the scattering length, but also 
on the energies of the colliding particles. The coupling parameter in one and two dimensions 
is therefore dependent on the chemical potential of the condensate. 

The energy dependent form of the many-body T-matrix in 2D found here can be used 
to obtain a self-consistent form for the 2D Gross-Pitaevskii equation. We have presented 
sample solutions and have shown that the importance of the nonlinear term is magnified in 
2D (as compared to the 3D case) due to the size of the coupling constant in two dimensions. 
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The Thomas-Fermi approximation is therefore valid at a much lower number of atoms than 
in the 3D case, approximately an order of magnitude lower in the case considered here. The 
critical frequency of vortex formation is also found to decrease with condensate occupation 
much faster in 2D than in 3D, and so vortices should be comparatively easier to form in 2D. 

Acknowledgments 

This research was supported by the Engineering and Physical Sciences Research Council of 
the United Kingdom, and by the European Union via the "Cold Quantum Gases" network. 
In addition, M.D.L. was supported by the Long Studentship from The Queen's College, 
Oxford, and S.A.M. would like to thank Trinity College, Oxford, for financial support. 



[1] E.B. Kolomeisky, T.J. Newman, J.P. Straley, and X. Qi, Phys. Rev. Lett. 85, 1146 (2000). 
[2] E.B. Kolomeisky and J.P. Straley, Phys. Rev. B 46, 11749 (1992). 

[3] A.I. Safonov, S.A. Vasilyev, I.S. Yasnikov, I.I. Lukashevich, and S. Jaakkola, Phys. Rev. Lett. 
81, 4545 (1998). 

[4] A. Gorlitz, J.M. Vogels, A.E. Leanhardt, C. Raman, T.L. Gustavson, J.R. Abo-Shaeer, 
A. P. Chikkatur, S. Gupta, S. Inouye, T. Rosenband, and W. Ketterle, Phys. Rev. Lett. 87, 
130402 (2001). 

[5] O. Zobay and B.M. Garraway, Phys. Rev. Lett. 86, 1195 (2001); H. Gauck, M. Hartl, 

D. Schneble, H. Schnitzler, T. Pfau, and J. Mlynek, Phys. Rev. Lett. 81, 5298 (1998); 

E. A. Hinds, M.G. Boshier, and I.G Hughes, Phys. Rev. Lett. 80, 645 (1998). 
[6] PC. Hohenberg, Phys. Rev. 158, 383 (1967). 

[7] J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 (1973). 

[8] V.N. Popov, Functional Integrals in Quantum Field Theory and Statistical Physics (Reidel, 
Dordrecht, 1983). 

[9] V. Bagnato and D. Kleppner, Phys. Rev. A 44,7439 (1991). 
[10] S.A. Morgan, J. Phys. B 33, 3847 (2000). 
[11] M. Bijlsma and H.T.C. Stoof, Phys. Rev. A 55, 498 (1997). 
[12] H. Shi and A. Griffin, Phys. Reps. 304, 1 (1998). 



21 



[13] D.S. Fisher and P.C. Hohenberg, Phys. Rev. B 37, 4936 (1988). 
[14] M. Schick, Phys. Rev. A 3, 1067 (1971). 
[15] B. Tanatar, Europhys. Lett. 51, 261 (2000). 
[16] G.S. Nunes, J. Phys. B 32, 4293 (1999). 

[17] D.S. Petrov, M. Holzmann, and G.V. Shlyapnikov, Phys. Rev. Lett. 84, 2551 (2000). 

[18] D.S. Petrov and G.V. Shlyapnikov, Phys. Rev. A 64, 012706 (2001). 

[19] S.A. Morgan, M.D. Lee, K. Burnett, Phys. Rev. A 65, 022706 (2002). 

[20] J.R. Taylor, Scattering Theory (Wiley, N.Y., 1972). 

[21] N.N. Bogoliubov, J. Phys. USSR 11, 23 (1947). 

[22] H.T.C. Stoof and M. Bijlsma, Phys. Rev. E 47, 939 (1993). 

[23] J.O. Anderson and H.T.C. Stoof, |cond-mat/0109143| . 

[24] H.T.C. Stoof, L.P.H. de Goey, W.M.H.M. Rovers, P.S.M. Kop Jansen, and B.J. Verhaar, Phys. 

Rev. A 38, 1248 (1988). 
[25] P. Bloom, Phys. Rev. B 12, 125 (1975). 
[26] A. Lenard, J. Math. Phys. 7, 1268 (1966). 
[27] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[28] M. Girardeau, J. Math. Phys. 1, 516 (1960). 

[29] N.P. Proukakis, S.A. Morgan, S. Choi, and K. Burnett, Phys. Rev. A 58, 2435 (1998); 
D.A.W. Hutchinson, K. Burnett, R.J. Dodd, S.A. Morgan, M. Rusch, E. Zaremba, 
N.P. Proukakis, M. Edwards, and C.W. Clark, J. Phys. B 33, 3825 (2000). 

[30] S.-H. Kim, C. Won, S.D. Oh, and W. Jhe, cond-mat/990487; S.-H. Kim, S.D. Oh, and W. Jhe, 
J. Korean Phys. Soc. 37, 665 (2000). 



22 



